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ABSTRACT 

We present resolved Herschel images of circumbinary debris disks in the a CrB (HD 
139006) and f3 Tri (HD13161) systems. By modelling their structure, we find that both 
disks are consistent with being aligned with the binary orbital planes. Though secular 
perturbations from the binary can bring the disk into alignment, in both cases the align- 
ment time at the distances at which the disk is resolved is greater than the stellar age, 
so we conclude that the coplanarity was primordial. Neither disk can be modelled as 
a narrow ring, requiring extended radial distributions. To satisfy both the Herschel and 
mid-IR images of the a CrB disk, we construct a model that extends from 1-300AU, 
whose radial profile is broadly consistent with a picture where planetesimal collisions 
are excited by secular perturbations from the binary. However, this model is also con- 
sistent with stirring by other mechanisms, such as the formation of Pluto- sized objects. 
The f3 Tri disk is modelled as a disk that extends from 50-400AU. A model with de- 
pleted (rather than empty) inner regions also reproduces the observations and is consis- 
tent with binary and other stirring mechanisms. As part of the modelling process, we 
find that the Herschel PACS beam varies by as much as 10% at 70/im and a few % 
at 100/im. The 70/im variation can therefore hinder image interpretation, particularly 
for poorly resolved objects. The number of systems in which circumbinary debris disk 
orientations have been compared with the binary plane is now four. More systems are 
needed, but a picture in which disks around very close binaries (a CrB, f3 Tri, and HD 
98800, with periods of a few weeks to a year) are aligned, and disks around wider bina- 
ries (99 Her, with a period of 50 years) are misaligned, may be emerging. This picture 
is qualitatively consistent with the expectation that the protoplanetary disks from which 
the debris emerged are more likely to be aligned if their binaries have shorter periods. 

Key words: circumstellar matter — stars: individual: f3 Trianguli, a Coronae Borealis, 
99 Herculis, HD 98800 



1 INTRODUCTION 

The Herschel Key Program DEBRIS (Dust Emission via a 
Bias free Reconnaissance in the Infrared/Submillimeter) has 
observed a large sample of nearby stars to discover and char- 
acterise extrasolar analogues to the Solar System's asteroid 
and Edge worth- Kuiper belts, collectively known as "debris 
disks." The 3.5m Hers chel mirror diamet er provides 6-7" res- 
olution at 70-100/xm (Pilbr attetalJ |2010), and as a conse- 
quence our survey has resolved many disks around stars in the 
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Solar neighbourhood for the first time (|Matthews et alJl20ldL 
Churcher et al ] |201ll : iKennedv et al1l2012L Wvatt et. al. 2012)! 1 ! 

Here we present resolved images of circumbinary disks in 
the a CrB and f3 Tri systems. These systems are interesting be- 
cause unlike most debris disk+binary systems, the binary orbits 
are well characterised. The combination of a known orbit and 
resolved disk means we can compare their relative inclination. 

Our observations of a disk around the binary 99 Her 



1 Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 
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( Kenn edy etaDl2012h wer e a step towa r d bu ilding on the bi- 
nary debris disk study of iTrilling et al J d2007h . Their Spitzer 
study found that debris disks are generally as common in bi- 
nary systems as in single systems, but are less likely to reside in 
systems with binary separation s in the 3-30AU range (see also 
iRodriguez & Zuckermani r2012). However, only some of their 
systems had detections at multiple wavelengths to constrain the 
disk location and none were reported as resolved, making the 
true dust location uncertain. Even in the case of dust detection 
at multiple wavelengths, the true dust location is unknown be- 
cause grains of different compositions and sizes can have the 
same temperature at different stellocentric distances. In addi- 
tion to uncertainty in the dust location, only the projected sky 
separation of the binary (not the binary semi-major axis) was 
generally known, adding further uncertainty. 

Systems with resolved disks and well characterised bi- 
nary orbits, such as 99 Her, a CrB, f3 Tri, and HD 98800 
dAndrews et al.ll201oh remove these ambiguities. The dust lo- 
cation and structure can be inferred within the context of the bi- 
nary orbit, leading to robust conclusions about whether the dust 
resides on stable orbits. One can also consider whether pertur- 
bations from the binary play an important role in setting the dust 
dynamics. For example, in the 99 Her system the disk position 
angle appears misaligned with the binary line of nodes. This 
misalignment is best explained with particles on polar orbits 
becau se these are stable for the stellar lifetime (iKennedv et al.l 
l2012h . Another question is whether binary perturbations can 
"stir" the disk particles by increasing their inclinations and ec- 
centricities, eventually resulting in high enough relative veloci- 
ties that collisions are destructive. Thi s process is analogous t o 
the planet stirring model proposed bv lMustill & WvatJ 12009). 
though may rely on vertical (inclination) stirring because com- 
panions of comparable mass induce lower forced eccentrici- 
ties than companions that are much less massive th an the star 
(IKennedv et alJl2012l : lMoriwaki & Nakagawall2004h . 

An additional link can be made to star and planet forma- 
tion. Young binary systems with small to medium (<100AU) 
separations are thought to form coplanar with their protoplan- 
etary disks, because the disk torque aligns the bina ry orbit on 
timescales short relative to the disk lifetime (e.g. iBate et all 
2000). Testing this prediction is difficult because it is difficult 
to ascertain disk and binary orien tations at the dista nces of the 
nearest star forming regions (e.g. lMonin et "aH l2007h . Because 
debris disks around older main- sequence binaries should retain 
the same orientation as the protoplanetary disks from which 
they emerged, these disks yield information on the outcome of 
star formation. The advantage is that compared to star-forming 
regions, these systems are much closer to Earth and hence larger 
on the sky and brighter, making disk and binary characterisation 
much easier. 

Though planet formation may be hindered to some de- 
gree by h igh collision velocities induced by binary per turba- 
tions (e.g. lMoriwaki & Nakagawall2004l : [Scholl et al.l2007h . the 
discovery of several circumbinary planets shows that planets 
do indeed form around binary stars (e.g. iDovle et alJl201ll : 
IWelsh etal.ll20l2h . Few such systems are known because bi- 
naries are generally avoided by ra dial velocity surveys (see 
Konacki 2005; Konacki et al. 2009). However, the very exis- 
tence of circumbinary debris disks provides evidence that planet 
formation around binaries can proceed to form at least 10- 
100km sized objects, w hich must exist to fee d the observed dust 
through collisions (e.g. Wyatt 2008; Krivovll2010h . Further, cir- 



cumbinary disks are relatively common dTrilling et al.ll2007h . 
suggesting that circumbinary planets may be no less unusual 
than their circumstellar equivalents. 

This paper is laid out as follows. We first consider the stel- 
lar and orbital properties, along with previous IR observations 
of the a CrB and (3 Tri systems (32]). We then show the Her- 
schel data ($3]) and simple models (^4]), and then interpret these 
within the context of the expected dynamics (33). We discuss 
the results and conclude in §3S]&[7] 



2 THE BINARY SYSTEMS 
2.1 a Coronae Borealis 

Among the "alpha" stars, a C rB (HD 139006 HIP 76267) is the 
only k nown eclipsing binary ds tebbinsll 1 9 1 4l : Tomkin & Popped 
119861) . Though some properties of this binary have been known 
for over a century fe.g. |jordanl [l910). only recently has the ori- 
entation of the orbit on the sky been constrained by the Hippar- 
cos mission (which als o provides a system distance of 23pc, 
IPerrvman & ESAlll997l : Ivan Leeuwenll2007h . The orbital ele- 
ments are given in Table T] The primary is an A O dwarf, or- 
bited by a G5 secondary (Tomkin & Popper 1 986h. The age of 
the system is about 350Myr ( Son g et al.l2QQll ; Rieke et al.l2005l : 
IVicanll2012h . 

T his system was discovered to have an IR excess with 
IRAS ( Aumann 1985), and the disk has subsequent ly been de- 
tected with the Spitzer M I PS0 and IRsB in struments dChen et all 
120061 : iRieke et all 120051 : ISuet all H0O6). The proximity and 
brightness mean this system has also been the target of high 
resolution mid-IR imaging campaigns, with resolution of the 
disk at 11 pm dlVIoerchen et al.|[201oh . They may have resolved 
the disk at 18/xm, but ambiguity arising from artefacts of the 
observing procedure meant that the results at this wavelength 
were unclear. The position angle of the extension in the 11/mi 
and first 18/xm images is roughly 350°, so consistent with the 
(rather uncertain) binary line of nodes. The images also suggest 
that the disk is edge on, so the inner part of the disk appears to 
be aligned with the binary orbital plane. 



2.2 p Trianguli 

The p Tri (HD 13161, HIP 10064) system was first recog- 
nised as a doub le lined spectroscopic binary over a century ago 
(Mitchell 1909). The orbital period was then found to be 37 
days, based on 12 ra dial velocity measurements. Further study 
has refined the orbit iStruve & Pogoll 19281 : lEbbighausenlll96(X) 
and modern instruments h ave since resolved th e pair allowing a 
visual orbit to be derived jHummel et all 1995b . The orbit is now 
well characterisedQ with the orbital element s given in Table [3 
The primary star has a A5IV spectral type dGrav et al]|2003h . 
and as the class IV indicates has probably reached the end of 
the main- sequence. The spectral type of the secondary is not 
known, but the mass suggests a mid F-t ype. There remain un - 
certainties in this system; as discussed bv lHummel et alj < fl995h . 



2 Multi-band Imaging Photometer for Spitzer (Werne r et al.1 12004 
IRieke et al.ll2004l) 

6 Infra-Red Spectrograph (Houck et al. 200|) 

4 The Washington Double Star Catalogue (WPS. iMason et al.ll200it) 
lists the orbit as "definitive" (grade 5) 
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Table 1. a C rB sys tem properties and la uncertainties 

iTomkin & Popped 1 19861 : IPerrvman & ESAl Il997h . The ascending 
node fl is measured anti-clockwise from North. The longitude of 
pericenter is measured anti-clockwise from the ascending node. The 
semi-major axis is calculated from the period, masses, and system 
distance. The luminosities are derived from our SED fitting to the 
AB pair (£]4.1K with the individual values found by scaling the 
Tomkin & Popper Il986h values of 74 and O.SLq. They derive a 
slightly greater system distance (26pc), which accounts for their larger 
luminosities. The uncertainty in La is set by the SED fi t ting, a nd 
for Lb we use the 25% uncertainty from lTomkin & Popper] dl986h to 
account for uncertainty in the flux ratio. 



Parameter 


Symbol (unit) 


Value 


Uncertainty 


Semi-major axis 


a (mas) 


8.66 


0.09 


Semi-major axis 


a(AU) 


0.2 


0.002 


Eccentricity 


e 


0.37 


0.01 


Inclination 


i(°) 


88.2 


0.1 


Ascending node 


Q(°) 


330 


20 


Longitude of pericenter 


"(°) 


311 


2 


Period 


P (days) 


17.3599 


0.0005 


Distance 


d(pc) 


23.0 


0.15 


M A 


(M ) 


2.58 


0.045 


M B 


(M ) 


0.92 


0.025 


La 


(£©) 


59.4 


1 


Lb 


(Lq) 


0.6 


0.15 



Table 2. (3 Tri system properties and la uncertainties (Pourbai3 l2000h . 
The ascending node Q is measured anti-clockwise from North. The 
longitude of pericenter is measured anti-clockwise from the ascending 
node. The total binary luminosity is derived from the SED in ^4.11 



Parameter 


Symbol (unit) 


Value 


Uncertainty 


Semi-major axis 


a (mas) 


8.03 


0.06 


Semi-major axis 


a(AU) 


0.312 


0.002 


Eccentricity 


e 


0.433 


0.004 


Inclination 


i(°) 


130 


0.52 


Ascending node 


no 


245.2 


0.67 


Longitude of pericenter 


w(°) 


118.1 


0.66 


Period 


P (days) 


31.3900 


0.0002 


Distance 


d(pc) 


38.9 


0.5 


M A 


(M ) 


3.5 


0.25 


M B 


(M ) 


1.4 


0.1 


Lab 


(£©) 


74 


1.4 



the A5 spectral type is at odds with the 3.5M mass, and the 
mass ratio is inconsistent wi th the luminosity rat i o. The system 
is simply classed as 'Old' bv lTrilling etaD(l2007h . lVicanl(l2012h 
derive an isochrone age of 730Myr, though this figure is sys- 
tematically uncertain both due to a relatively small luminosity 
ratio (so the age may be overestimated) and because (3 Tri lies 
outside the range encompassed by some isochrone models. The 
exact age is not particularly important for our purposes here, so 
we assume it is 730Myr. 

Like a CrB, t his system was first discove red to have an IR 
excess with IRAS dSadakane & Nishidall 19861) . It was observed 
with Spitzer MIPS as p art of a sample of main- sequence bina- 
ries (iTrilling et al.ll2007h . and also observed with IRS in the two 
longer wavelength modules (14-38/im). 
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Table 3. Herschel observations of a CrB and (5 Tri. 



Target 


Obsld 


Date 


Instrument 


Duration (s) 


a: CrB 


1342213794 


7 Feb 2011 


PACS 100/160 


445 


a CrB 


1342213795 


7 Feb 2011 


PACS 100/160 


445 


(3 Tri 




4 Tnlv 901 1 


PArN 100/160 


445 


/3Tri 


1342223651 


4 July 2011 


PACS 100/160 


445 


a: CrB 


1342223846 


9 July 2011 


PACS 70/160 


445 


a: CrB 


1342223847 


9 July 2011 


PACS 70/160 


445 


/3Tri 


1342237390 


12 Jan 2012 


PACS 70/160 


445 


(3 Tri 


1342237391 


12 Jan 2012 


PACS 70/160 


445 


/3Tri 


1342237504 


14 Jan 2012 


SPIRE 250/350/500 


721 



3 OBSERVATIONS 
3.1 Herschel 

Hersche l Photodetector and Array Camera & Spectrometer 
(PACS, Podit sch et all hoiO) data were taken for both a 
CrB and /3 Tri at 100 and 160/mi during routine DEBRIS ob- 
servations. Subse quently, a Spectral and Photometric Imaging 
Receiver (SPIRE. [Griffin et"aUl2010l) observation of /3 Tri was 
triggered by the large PACS excess indicating a likely sub-mm 
detection. For both targets we also obtained 70/im PACS images 
to better resolve the disks. Because every PACS observation in- 
cludes the 160/im band, we have two images at this wavelength 
for each target. All observations were taken in the standard scan- 
map modes for our survey; mini scan-maps for PACS data and 
small maps for SPIRE. Data were reduced using a near-standard 
pipeline with the Hersc hel Int eractive Processing Environment 
(HIPE Version 7.0. lot3l2oTob . We decrease the noise slightly 
by including some "turn-around" data taken as the telescope is 
accelerating and decelerating at the start and end of each scan 
leg. 

Figure Q] shows the Herschel PACS images of a CrB. The 
disk is not obviously resolved, so the inset shows the residuals 
from a point- spread function (PSF) fit to the 100/xm data (an ob- 
servation of Herschel calibrator 7 Dra generated using the same 
data reduction method was used as a PSF, see below). The resid- 
ual emission in the inset is symmetric about the stellar position 
along a position angle of about 350°, showing that the disk is 
most likely resolved with a near edge-on geometry. The emis- 
sion extends to around 10" (^230AU) from either side of the 
stellar position, suggesting a disk diameter of several hundred 
AU. The disk does not appear to be resolved along the minor 
axis, so is consistent with the disk being edge on and aligned 
with the (eclipsing) binary plane. A PSF fit to the 70/im image 
shows similar results, while the 160/im image is unresolved. 

Because the a CrB disk is not well resolved, it is possi- 
ble that this extension is affected by variation in the PACS PSF 
(or beam). Further, a key difference between our method and 
that used for beam characterisation by the PACS team is that 
because they are much fainter than the calibration stars our data 
cannot be re-centered on a frame-by-frame basis, resulting in 
a slightly larger beam size. Thus, comparison of the Gaussian 
size in our images with those reported by the PACS team could 
lead to the conclusion that the disk is slightly more extended 
than it really is. To characterise the beam variation we have ob- 
tained two observations each of the five PACS calibration stars 
at both 70 and 100/xm. These observations were nearly all taken 
using the same mini scan-map mode used by DEBRIS. These 
data are reduced in the same way as all DEBRIS data, so allow 
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Figure 1. Resolved PACS images of a CrB at 70 (left) and 100/im (right). North is up and East is left, and the arrow shows the spacecraft "Z" 
direction. The colour scale is in mJy/square arcsecond. The overlaid contours show the corresponding 160/im images in 5 linear steps from 5-13cr. 
Similar contours correspond to 30-200(7 for the 70/im image, and 10-75(7 for the 100/im image. The hatched circles show the average PACS beam 
FWHM of 5775 (70/im) and 6 / /87 (100/im, see text). The inset (25" square) shows residuals of the 100/mi a CrB image after PSF fitting, the lobes 
show that the disk is clearly resolved. 



a realistic comparison. The expectation is that each calibration 
PSF will be slightly different, with the differences across all ten 
being representative of the uncertainty in the PSF specific to a 
given science observation. 

The 70/im calibration observations are shown in Figure |3 
where each image shows the residuals when the average of all 
ten observations is subtracted (after the peak is scaled to the 
same as the observation). All images have been rotated so that 
spacecraft "Z" is up (spacecraft "Y" is to the right). Each panel 
shows the observation number (ObsID), star name, and Gaus- 
sian major and minor FWHM. The ellipse shows the Gaussian 
FWHM and the line shows the major axis orientation. Because 
the expectation is that each observation is independent, the pan- 
els are ordered chronologically (i.e. by ObsID). It is clear that 
the beam width varies at about the 10% level along the minor 
axis, between 5728 for a Boo and 5 / /79 for a Cet. The vari- 
ation along the major axis is much smaller at about 2%. The 
beam variation appears to be systematically different for two 
stars; both the a Boo and a Tau images show similarly small 
minor FWHM despite being observed about six months apart. 
This variation cannot be attributed to circumstellar material be- 
cause none of the calibration stars have IR excesses. It may be 
difficult to confirm this possible systematic effect, for example 
as a function of Solar elongation or declination, as there are only 
five PACS calibration stars. 

The major and minor FWHM of the average 70/im PACS 
beam are 5'/ 91 and 5759, so similar to, but slightly larger than, 
the values of 5776 and 5746 found by the PACS team (for a 
scan speed of 20" /sec)0 Based on the variation seen in Fisure 
13 the uncertainties on these values are about 0705 in the ma- 
jor axis (approximately along the spacecraft "Y" direction) and 



5 See PACS Observer's Manual 



072 in the minor axis (approximately along the spacecraft "Z" 
direction). 

A similar analysis for the 100/im beam is shown in Figure 
[3] The panels are again in chronological order, so the order of 
the stars is not the same as at 70/im. The variation is smaller 
at this wavelength (2-4%), but some systematic difference for a 
Tau and a Boo is still seen. The average major and minor beam 
FWHM are 6795 and 6778 (compared to 6 / /89 and 6 / /69 found 
by the PACS team), with uncertainties of about 0" 1 in both axes. 
The lower left panel (a Cet, ObsID: 1342203033) shows rather 
different residuals compared to all others, arising because the 
observation comprises only one scan direction. This difference 
shows that aside from variation, the beam is also a function of 
the observing strategy, and that the strategy should be the same 
when using calibration stars to analyse science data (e.g. for 
PSF fitting and image modelling). 

Having quantified the beam characteristics specific to our 
data, we now compare them to the a CrB observations. A Gaus- 
sian fit to the star- subtracted 70/im image finds a position angle 
of 348 =b 3°, which is consistent with the ascending node of 
330 db 20° for the binary. The major and minor FWHM are 
7" 3 ± 07l and 67 1 ± 07l. Figure shows that a major axis 
of the beam is typically perpendicular to the spacecraft Z axis, 
so the disk position angle is between the major and minor axes, 
where the beam has a FWHM of about 5.75 ± 07l. Therefore, 
the disk is clearly resolved along the direction of the PA, but at 
only about 2a significance in the perpendicular direction. Sim- 
ple deconvolution suggests a disk diameter of about 100AU. 
This size is much smaller than suggested by the residuals in 
the inset in Figure [T] indicating that the bulk of the emission 
is poorly resolved. This structure is perhaps consistent with the 
mid-IR imaging if a CrB hosts warm and cold components that 
both contribute to the Herschel fluxes but only the cold com- 
ponent is resolved with Herschel. Though the minor Gaussian 
FWHM is slightly larger than expected, variation in the 70/im 
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Figure 2. PACS 70/im beam comparison. Each panel shows the beam after the average beam has been subtracted, with labels noting ObsID, star 
name, and Gaussian major and minor FWHM. The ellipse shows the FWHM, with the line along the major axis. The colour scale shows the fractional 
difference (in %) relative to the peak. In these images the spacecraft "Z" direction is up and "Y" is to the right, and the two scan directions that make a 
mini scan map are ±20° from horizontal. 
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Figure 3. PACS 100/im beam comparison. The layout is as described in Figured 



beam means that the disk is consistent with being edge-on, and 
therefore aligned with the (eclipsing) binary orbital plane. This 
analysis also shows that the outer component resolved by Her- 
schel is consistent with being aligned with the i nner component 
resolved in the mid-IR dMoerchen et alJl201(t) . At 100/zm the 
major and minor FWHM are 8"2 ± 0" 16 and 6 x /9 ± Of! 13, yield- 
ing a size of about 100AU and a position angle (PA) of 347 ±4° . 
Therefore, at this wavelength the disk is resolved along the disk 
PA, but not in the perpendicular direction, so is again consistent 
with being edge-on, and aligned with the binary plane. We re- 
turn to the issue of beam variation when creating resolved mod- 
els of the a CrB disk in fl4~2l 

The f3 Tri images are shown in Figure HI Compared to 
the beam size, the disk around (3 Tri is clearly resolved at 70 
and 100/im. The disk is also resolved at 160/xm, but not at the 
SPIRE wavelengths of 250-500/xm. A Gaussian fit to the star- 
subtracted image at 70/im finds major and minor FWHM of 
9'/ 1 ± 0'/ 13 and fl 5 ± 07 1 at a PA of 66 ± 3° . Simple Gaussian 



deconvolution suggest a disk size of about 270AU, and an incli- 
nation of about 46 ± 3° . At 100/xm the FWHM are 1074 ± 07 13 
and 8 x /6 =b O'll at a PA of 68 =b 3°, suggesting a size of about 
300AU and inclination of 48 ± 3°. At 160/xm the FWHM are 
14 ; /4 ± 0712 and 1273 ± 07l at a PA of 68 ± 2° (the 160/xm 
beam is elongated and about 1077 by 127 1). The 160/xm image 
is therefore poorly resolved, but suggests an approximate size 
of 300-370AU. Taken together, these angles mean that in this 
system the disk is again consistent with being aligned with the 
binary plane. The increasing disk size with wavelength suggests 
that the disk may be extended. 

The flux density in the PACS images for both a CrB and /3 
Tri is measured using apertures and the SPIRE fluxes with PSF 
fitting. The measurements are given in Table [4] The uncertain- 
ties include both statistical and s ystematic uncertaintie s due to 
repeatability and calibration (see iKennedv et al.ll20li for fur- 
ther comments on calibration). 
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Figure 4. Resolved PACS images of (3 Tri at 70 (left) and 100/im (right). North is up and East is left. The colour scale is in mJy/square arcsecond. 
The overlaid contours show the corresponding 160/im images in 5 linear steps from 5-30cr. Similar contours correspond to 5-75cr for the 70/xm image, 
and 5-70(7 for the 100/xm image. The hatched circles show the average PACS beam FWHM of 5'!l5 and 6 / /87 (see text). The inset (25" square) shows 
residuals of the 100/xm (3 Tri image after PSF fitting, the elliptic residual emission shows that the disk is clearly resolved. 



Table 4. Herschel photometry and la uncertainties of a CrB and (3 Tri. 
The 160/im measurements from each observation are given separately. 



Hauck & Mermi lliodll 19971: IPerrvman & ESA|[l997l: iHog et al I 



Target 


Band 


Flux (mJy) 


Uncertainty (mJy) 


Method 


a CrB 


PACS70 


515 


26 


15" aperture 


a CrB 


PACS 160 


69 


5 


PSF fit 


a CrB 


PACS 100 


235 


12 


15" aperture 


a CrB 


PACS 160 


66 


5 


PSF fit 


(3 Tri 


PACS70 


641 


32 


20" aperture 


(3 Tri 


PACS 160 


264 


17 


25" aperture 


(3 Tri 


PACS 100 


481 


23 


20" aperture 


(3 Tri 


PACS 160 


264 


16 


25" aperture 


(3 Tri 


SPIRE250 


87 


7 


PSF fit 


(3 Tri 


SPIRE350 


35 


5 


PSF fit 


/3Tri 


SPIRE500 


5 


5 


PSF fit 



3.2 Spitzer 

As noted above, both systems were observed by Spitzer with the 
MIPS and IRS instruments. The MIPS photometry is taken from 
Phillip] ([201 1), and the IRS spectra retrieved from the Cor- 
nell Atlas of Spit zer/Infrared Spectrograph Sources (CASSIS, 
iLebouteiller et al.ll201ll) . The PACS and MIPS 70/xm measure- 
ments are consistent (460 =b 47mJy for a CrB and 677 =b 62mJy 
for f3 Tri). Though we instead use the PACS images, f3 Tri is 
clearly resolved at 70/xm with Spitzer. 



4 MODELS 
4.1 SEDs 

We first model the stellar and disk photometry to derive 
approximate disk temperatures and fractional luminosi- 
ties. This photometry is c ol lected from numero u s cat- 
alogues dMorel & Magnenatl Il978l : Moshir etjd] 1 19931 : 



20001: ICutri et al.ll2003l: lMermilliodll2006l : llshihara et alfeoiOl : 
Phillips 201 ll: ILebouteiller et al.l201ll) . Photometry and colours 
at wavelengths up to 9/xm are used to mo del the stellar pho- 
tosphere, using the PHOENIX Gaia grid terott & H auschildt 
120051) . The best fit model is found by least squares min- 
imisation. At wavelengths longer than 9/xm we model the 
excess emission above the photosphere with a simple modified 
blackbody function; at wavelengths beyond Ao the blackbody 
is multiplied by (A/Ao) _/3 . 

The spectral energy distributions (SEDs) are shown in Fig- 
ure \5\ The best fitting stellar model for a CrB has T e ff = 
9280 ± 100K, R* = 3.06 ± O.O3i? , and L* = 60 ± LL , 
and for f3 Tri has T eff = 8000 ± 100K, R* = 4.6 ± O.O5i? , 
and L* = 74 ± 2L . Because no photometry resolves either 
binary, these parameters simply describe what the total stellar 
emission spectra look like, and are not physical. The a CrB disk 
model has temperature Tdisk = 124K, fractional luminosity 
/ = 1.7 x 10 -5 , corresponding to a total grain surface area 
cr tot — 0.34AU 2 for blackbody grains, and Ao = 75/xm and 
(3 = 1.6. The f3 Tri model has temperature Tdisk = 84K, 
/ = 3.0 x 10~ 5 , total grain surface area cr to t = 3.3AU 2 for 
blackbody grains, and Ao = 152/xm and f3 = 1.25. 

Despite the sugge stion of warm dust by mid-IR imaging 
dMoerchen et al.ll2010h . the a CrB disk SED, most notably the 
IRS spectrum, does not require dust at multiple temperatures. 
The inferred disk radius assuming blackbody properties for the 
single component fit to the SED is 40AU, well beyond the few 
AU scale of the mid-IR emission. The SED does not strongly 
preclude the existence of warm dust however, because it is lim- 
ited by the ~2% level that is achievable with typical photomet- 
ric calibration, and a similar uncertainty for the IRS spectrum. 
Therefore, at 10-20/xm an excess of <100mJy could be present 
yet not confidently detected. How this limit converts to a limit 
on the fractional luminosity depends on the dust temperature 
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Figure 5. SEDs for a CrB (left) and (5 Tri (right). Photometry is shown as black dots or black triangles for upper limits. Disk (i.e. photosphere- 
subtracted) fluxes and upper limits are shown as grey dots, open triangles, and small green dots for IRS. The stellar spectrum is shown as a blue line 
and the blackbody disk model as a red line, with the total shown as a black line. 



(e.g. lWvattl2008b . In the case of 300K dust the limit on the frac- 
tional luminosity is about 10 -5 , so only about a factor of two 
lower than derived from the 124K disk SED above. 

For f3 Tri the blackbody radius is larger at around 95 AU. 
This size is somewhat smaller than the 150AU derived from the 
Gaussian fitting, but as noted above the disk could be extended. 
Aside from being marginally resolved with MIPS, there are no 
constraints on the structure from other observations. 



4.2 Images 

The main goal of modelling the resolved images is to derive 
the orientation of the disk in space, for comparison with the 
orientation of the binary orbital plane. Though alignment is al- 
ready suggested by Gaussian fitting in both cases, the PACS 
beam is not azimuthally symmetric so resolved modelling pro- 
vides a worthwhile check. To achieve this goal, we must find a 
spatial structure that reproduces the images satisfactorily. The 
model s are gener ated using the metho d described in lWvatt et al.1 
d2012h (see also iKennedv et al.ll2012l) . Basically, a high reso- 
lution model of the structure is generated and multiplied by 
the grain emission properties at each wavelength. The struc- 
ture is the spatial distribution of dust, implemented as the cross- 
sectional area as a function of radial distance from the star. The 
distribution is disk-like with a small range of particle inclina- 
tions and the dust's face-on optical depth (r) is parameterised as 
a function of radius by one or more power laws. The grain emis- 
sion properties are simply described with a modified blackbody 
whose temperature decreases with radial stellocentric distance 
as a power-law. The high resolution models are then created 
with some spatial orientation and convolved with the instrument 
beam for comparison with the observed images. The large num- 
ber of model parameters, as well as different possible configura- 
tions (e.g. extended disks vs. multiple rings) means that we do 
not undertake a grid search of possible parameter spaces. The 
final models are found by a combination of by-eye fitting and 
least- squares minimisation. 

Our general approach to modelling is to use the simplest 
model possible to explain the data in hand, increasing the model 



complexity as required. Here, we find that the simplest model, 
a narrow ring, is not sufficient to reproduce the PACS images 
for either system. As a second step we use extended dust dis- 
tributions, though two discrete narrow rings is also a possibility 
(e.g. Wyatt et al. 2012), as is the inclusion of an unresolved 
component. Because we do not use a grid approach, estimating 
uncertainties on model parameters, many of which are highly 
degenerate, is made more difficult, However, we find that mul- 
tiple different configurations can reasonably reproduce the data, 
which gives an idea of how well the data constrain the disk 
structure given the relatively poor resolution compared to the 
disk sizes. Similarly, we find no requirement for the temper- 
ature distribution to deviate from a simple blackbody relation 
(T = 278.3Ll /4 r _0 5 K, where L* is the binary luminosity in 
Solar units, and r is the disk radius in AU). We do not claim 
this is the true temperature distribution, given that grain temper- 
atures are expected to be somewhat different (hotter) than the 
blackbody relation suggests, but that there is insufficient infor- 
mation to quantify this difference. 



4.2.1 a CrB 

We first model a CrB, for which about half to two thirds of 
the Herschel emission is unresolved, with the remainder more 
distant from the stellar position and resolved (see Fig. [6]). This 
large level of unresolved flux means that different models of the 
system can satisfy the observations, with the main constraint 
being that mid-IR emission must be detected only within a few 
AU of the star. Because this detection could be the Wein side 
of the emission, it does not preclude the existence of material 
just beyond a few AU. That is, the mid-IR emission does not 
indicate that no dust exists where it was not detected, because it 
could be cool enough to evade detection. 

The size of th e 11 /mi images (described in ^2.11 
Moerchen et al .2010) compared to the reference PSF suggests a 
disk radius of about 2.3AU. Given that this warm emission can- 
not be resolved by Herschel, there must either be (at least) two 
separate disk components (i.e. inner and outer components), or 
a continuous disk that extends from a few to a few hundred AU, 
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Figure 6. Residuals after the resolved model of a CrB is subtracted from the observations at 70/im (left) 100/im (middle) and 160/im (right). The 
colour scale is in units of the pixel to pixel uncertainty, and contours are shown at ±3, 4 and 5a. The white contour shows the 3 a contour from the 
observed images. 



whose structure appears different at different wavelengths due 
to the changing temperature with radial distance. 

As an example, the Herschel data and SED can be mod- 
elled with an outer disk that extends from 45-1 95 AU with a 
decreasing face-on optical depth r = 2.1 x 10 -5 r -2-9 , and 
an unresolved inner component at a temperature of 120K. The 
outer component has a position angle of 345° and is edge-on, 
so consistent with being aligned with the binary orbital plane. 
The interpretation of this model is that the unresolved emission 
is distributed such that some dust is close enough to the star to 
provide the mid-IR emission and structure, while most of it lies 
distant enough from the star so that overall it appears to have 
a temperature of 120K. However, the distance inferred for a 
blackbody at 120K around a CrB is 40AU, or l'/8, which is rel- 
atively large given the 6-7" FWHM PACS beam at 70-100/xm, 
and should be detectable in the images. 

We therefore use a more complex two power-law model, 
in which the disk extends continuously from very near the star 
to several hundred AU. The surface density increases from the 
inner edge to some break radius, beyond which it decreases 
to the outer edge. This model has the advantage that it is ex- 
tended near the binary, so may be able to explain the mid- 
IR emission without invoking a separate warm component. 
Such a surface density structure may be expected, because 
debris disk deca y (and stirring in some models) is an inside 
out process (e.g.lKenvon & Bromlev 20 041: IWyatt et al .Il2007al : 



iKennedv & Wvattll2010l : 



Mustill & Wyatt 2009). 



Figure [6] shows the results of implementing this model, 
where the dust extends from 1AU to 300 AU, with the break 
radius at 50AU. The inner edge of 1AU corresponds to the in- 
ner li mit of stability for parti cles orbiting the a CrB binary 
(e.g. Holman & Wiegert 1999). The optical depth of the inner 
component is r = 6.4 x 10 _8 r 17 . For the outer component 
r = 3.3 x 10 -3 r -1-6 . The disk temperature has a blackbody 
dependence, with Tdisk = 775r _0-5 K. As can be seen from 
the SED (Fig. [5}, some modification of the blackbody function 
is needed, with Ao = 69/xm and /3 = 1.3 for the inner com- 
ponent and Ao = 82/xm and f3 = 2.2 for the outer compo- 
nent. The poor resolution relative to the disk size means that 
these parameters are poorly constrained with strong degenera- 
cies. There is a factor of eight drop in surface density at the 
break radius; attempts to keep the dust distribution continuous 
within the two power-law model were unsuccessful because the 



outer component is fairly extended but with low surface bright- 
ness. Joining the outer component smoothly to the inner com- 
ponent results in too much flux just beyond the break radius. 
Such a drop is an expected featu re of some disk evolution mod- 
els (e.g. Kennedy & Wyatt 2010), a point we return to in ^6] The 
extent of the outer component is not very well constrained be- 
cause the surface brightness decreases strongly with radius. The 
disk model has a position angle of 345° so is consistent with 
the binary line of nodes. The disk is edge on, so again consis- 
tent with being coplanar with the binary orbit. The images show 
significant residuals for changes of 10-20° in both inclination 
and PA, giving an indication of their uncertainties. 

The two models considered above show similar residual 
structure when the model is subtracted from the observations. 
The only sizeable (>3cr) residual related to the star+disk is 
located on one side of the stellar position at 70/xm. With the 
simple models considered here, similar residuals remain for a 
range of PSFs (i.e. those shown in Fig. [2}, even if the open- 
ing angle of the inner disk component is allowed to vary (as 
might be expected if this component was misaligned with the 
outer component, see ^5] below). Our inability to construct a 
model that accounts for these residuals leads us to conclude that 
they are likely due to the beam variations discussed in ^3] Al- 
ternatively, there may be additional structure unaccounted for 
by our model, which would have a scale of a few arcseconds 
to leave such small residuals. While the mid-IR images suggest 
that our coplanar configuration is representative, these observa- 
tions should be repeated, specifically at 18/xm, as a further test 
of disk coplanarity and structure near the star. 

To compare the e xtended model with the resolved structure 
seen in the mid-IR by iMoerchen et al.l (120101) . Figure [7] shows 
synthetic images at 11 and 18/xm, where we have convolved 
the star+disk models with appropriately sized Gaussian PSFs, 
added noise, and then subtracted t he PSFs scaled to the image 
peak (i.e. the same method used bv lMoerchen et al.l2010[ to de- 
tect extension). The background noise level is calculated from 
the photometric uncertainties of 2 and 15mJy at 11 and 18/xm 
respectively, assuming aperture radii of 2 / /0. In these images we 
have multiplied our model by a factor of two at 11/xm, and five 
at 18/xm to make them appear similar to the mid-IR observa- 
tions, and show how the observed mid-IR structure could arise 
from an extended disk. A factor of two at 11 /im results in a 
disk flux of 65mJy, so is reasonable given the uncertainty in the 



© 0000 RAS, MNRAS 000, 000-000 



Coplanar Circumbinary Debris Disks 9 




2 10-1-2 210-1 

a cos 8 offset (") a cos 8 offset (") 



Figur e 7. Simulated 11 /im (left) and 18/xm (right) peak- subtracted mid-IR observations of the a CrB model for comparison with Moerchen et al. 
d2010h . The pixel scale is 0.1'Vpixel, and the scales are in mJy/pixel. Contours are drawn at 3, 6, and 9 times the background noise level. Our model 
has been multiplied by a factor of two at 1 1/xm, and five at 18/xm (see text). 



disk flux at this wavelength, which is limited by the calibration 
of the photosphere and IRS spectrum. The factor five at 18/xm 
results in a disk flux of lOOOmJy, which is less reasonable and 
would make the disk flux inconsistent with the IRS and AKARI 
measurements. Our disk model is therefore fainter at 18/im than 
could be detected with the mid-IR imaging. That the 18/xm ob- 
servations suffer from artefacts due to the observing procedure 
is the likely explanation, providing further motivation for re- 
peating these observations. 



4.2.2 p Tri 

The simplest (3 Tri model is shown in the top row of Figure 
[8] We found that the images could not be modelled as a sin- 
gle narrow ring, so first allowed the inner and outer disk edges 
to vary with a power-law surface density profile in between. 
This model extends from 50-400AU with an optical depth of 
r — 6.7 x 10 -5 r -0-9 . The temperature dependence is again 
that for a blackbody (T dis k = 819r~ 05 K), with A = 130/xm 
and f3 = 1.5. The disk outer edge is not well constrained by the 
PACS images, but is constrained to some degree by not being 
resolved with SPIRE at 250/xm. The position angle and inclina- 
tion are 67° and 46° respectively. Variation of the inclination by 
10° leaves significant residuals compared to the data, as does a 
similar change in the position angle. The position angle of the 
binary plane is 65° so is easily consistent with that found for the 
disk. The inclination of the binary plane is 40°, so not signifi- 
cantly different to that found for the disk. While the inclination 
is fairly well constrained the disk opening angle is not, with val- 
ues between 0-50° giving similarly good fits (though lower disk 
opening angles are preferred). As we discuss in detail below in 
§5\ this ambiguity can be resolved by considering circumbinary 
particle dynamics. 

As with a CrB, we also tried a continuous surface den- 
sity profile by adding an inner component to the first model. 
This continuous model extends from the innermost stable or- 



bit at about 4AU lHolman & Wiegerj dl999h out to 500AU with 
the surface density peaking at 100AU and is equally successful 
at reproducing the Herschel images (bottom row of Fig Hi. In 
this model we restricted the inner component to join smoothly 
to the outer component; without this restriction the inner com- 
ponent can be arbitrarily small and the model no different to 
the extended one above. The inner component has an optical 
depth ofr = 2 x 10 _8 r 17 , and the outer component decreases 
with an index of -1.7. The temperature is for a blackbody as 
above, with Ao = 85/xm and f3 — 1 for the inner component 
and Ao = 137/xm and ft = 1.5 for the outer component. As 
with a CrB, these parameters are poorly constrained and highly 
degenerate. 

While these two models show that the radial and vertical 
disk structure is poorly constrained, both suggest that the disk 
optical depth decreases beyond a maximum that lies around 60- 
100AU. Interior regions may be depleted as in the continuous 
model, or completely empty. Despite these uncertainties, our 
main conclusion that the disk is consistent with being aligned 
with the binary orbital plane is robust. 



5 DYNAMICS 

The main result of the resolved modelling is that the disks are 
both consistent with being aligned with the orbital planes of the 
host binaries. Whether this alignment simply reflects the initial 
debris disk+binary configuration at the end of the protoplane- 
tary disk phase, or subsequent evolution requires some study of 
the expected dynamics. 

Disk particles exterior to a pair of orbiting bodies (in our 
case a binary star system) are subject to secular perturbations 
that modify the particle orbits on long timescales. How the par- 
ticle orbits are modified depends on the properties of the bi- 
nary, particularly the mass ratio and binary eccentricity, but 
also depends on the particle o r bits themselves (jVerrier & Evansl 
120091 : iFarago & LaskaJbOld : iDoolin & BlundelllbOllh . Ifthe 
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Figure 8. Residuals after the resolved model of (3 Tri is subtracted from the observations at 70/xm (left) 100/xm (middle) and 160/xm (right). The top 
row shows results from the extended (50-400 AU) model, and the bottom row shows results from the continuous (1-400AU) model. The colour scale is 
in units of the pixel to pixel uncertainty, and contours are shown at ±3 a. The white contour shows the 3a contour from the observed images. 



stars in the binary have similar masses, the effect on par- 
ticle eccentricities is smalle r than for larger m ass ratios 
(iMoriwaki & Nakagawal 120041 : iKennedv et alTl2012h . In con- 
trast, the changes in particle inclinations and nodes can be sig- 
nificant. If the angle between the particle and binary orbital 
planes is small (i.e. particles have low inclinations with respect 
to the binary plane, we quantify "small" and "large" inclinations 
below), particle inclinations oscillate about the binary plane. 
Therefore, while a disk of such particles will have an opening 
angle double the initial misalignment, it will also appear to be 
aligned with the binary plane. The best example is the warp in 
the f3 Pictoris disk, which was proposed to be d ue to perturba- 
tions by a planet misaligned with the disk plane dMouillet et al.l 



1 1997; IXugereau et all 
(lLagrange et al.ll2009l 



2001) that was subsequently discovered 
201(3). For large differences between the 



particle and binary planes (i.e. high inclinations with respect to 
the binary plane), particle inclinations are coupled to the evo- 
lution of their line of nodes, and oscillate around a polar orbit. 
In the context of a debris disk these particles evolve into vari- 
ous structures depend ing on the degree of initial misalignment 
(Kenn edy etaDl2012h . Within the high-inclination family of or- 
bits, only disks that are initially misaligned by about 90° and 
have their line of nodes perpendicular to the binary pericenter 
will not be strongly perturbed and continue to appear disk-like. 

The dividing line between "large" and "small" relative in- 
clinations, and thu s the different families of disk structures, 
was quantified by iFarago & Laskarl d201Qh . For a CrB with 
e = 0.37 the critical angle is 48°, for f3 Tri with e = 0.433 
the critical inclination is 43° (for an ascending node of ±90°). 
Thus, if the disk inclination relative to the binary plane was ini- 



tially larger than 48 and 43° respectively for these systems, it 
could not become aligned with the binary plane. 

While particles with initial inclinations lower than the crit- 
ical angle can become aligned, the finite secular precession time 
means that they will only be aligned if they have been perturbed 
over a suf ficiently long period. The timescale for alignment is 
given bv lFarago & Laskarj (120101) . and is shown for an arbitrar- 
ily small initial disk-binary plane inclination, and for initial mis- 
alignments of 20 and 40° for a CrB and f3 Tri in Figure|9] Parti- 
cles with semi-major axes that lie to the left of where the curves 
intersect a given system age (i.e. less than r a ii gn ) have com- 
pleted at least one cycle of secular evolution. The hatched re- 
gion indicates the region where PACS cannot resolve the disk at 
70/im, so shows where the disk structure is poorly constrained 
by Herschel observations. 

The a CrB system is about 350Myr old, so the inclinations 
of particles within about 80AU will be symmetric about the bi- 
nary orbital plane. Beyond this distance however, particles have 
not yet had time to undergo a complete cycle of secular preces- 
sion and will retain their original orbital inclinations. The tran- 
sition between these regions is n ot sharp due to the fin ite secu- 
lar precession time (see Fig. 1 of Augereauet Figure 
[9] also shows that the alignment distance is about the same as 
the inner extent of the PACS resolution at 70/xm. Thus, the disk 
structure that is constrained by the PACS images lies beyond 
Malign, and the alignment with the binary plane cannot arise due 
to secular perturbations. This conclusion is strengthened by the 
finding that the orientation as seen by Herschel is consistent 
with the inner disk as seen in the mid-IR. 

In the f3 Tri disk, particles beyond r a ii gn ~ 140AU are 
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Figure 9. Secular precession times for a CrB and (3 Tri (solid lines, for misalignment angles of 0, 20, and 40° from bottom to top). The dot dashed 
line shows the estimated stellar age, and the dotted lines disk radii (the break and outer radii for the continuous a CrB model, and the inner and outer 
radii for the extended (3 Tri model). The alignment radius is where the secular precession time equals the stellar age for low misalignment angles. The 
hatched region lies inside the PACS beam half- width half-maximum at 70/im, so structure in this region is poorly constrained by PACS observations. 



too distant to be significantly affected by secular perturbations 
from the binary, assuming an age of 730Myr. Because we find 
that the disk extends well beyond this distance, the alignment 
of the binary orbital and disk planes are again most likely pri- 
mordial. Though we could not constrain the disk opening angle 
well, the assumption that it is small is reasonable based on i) ob- 
servations that edge-on debris disks typically have small open- 
ing angles (e.g. lKrist etaDl2005l : iGolimowski et alJl2006h . and 
ii) that secular perturbations could not have increased the scale 
height at the distances resolved by the Herschel observations. 



6 DISCUSSION 

By considering the resolved models and expected dynamics, 
we have shown that the circumbinary debris disks in both the 
a CrB and (3 Tri systems probably formed coplanar with their 
parent binaries. A corollary is that because the disks were pri- 
mordially aligned, there should be little vertical structure inside 
the alignment radius. Both disks were successfully modelled as 
disks with a single plane of symmetry, but are too poorly re- 
solved with Herschel to strongly verify this statement. However, 
for a CrB the low-level contours from the 11 and 18/xm mid-IR 
imaging suggest that the inner disk has a similar position an- 
gle to the binary and outer disk. Only one of the 18/xm images 
shows the same position angle, though the fact that the other is 
actually narrower than the PSF in the s ame direction makes thi s 
extension questionable, as discussed bv lMoerchen et al.l d2010h . 

While there appears to be no need or evidence for vertical 
disk structure induced by perturbations to disk particles' incli- 
nations, an indirect signature could exist from inclination and 
eccentricity variations imposed on disk particles. These varia- 
tions "stir" the disk, where we are defining "stirring" to be any 
mechanism that increases relative velocities between particles 
sufficiently that collisions become catastrophically destructive, 
thus initiating a collisional cascade. In stirred regions mass at 
the small end of the collisional cascade is removed by radiation 
forces, thus depleting the disk. In unstirred regions the mass re- 
mains constant. The size distribution in stirred regions has many 



more small grains for a given mass in large objects, so is more 
easily detected due to the larger emitting surface area0 

Though the structure of the inner disk is not well con- 
strained, we have shown that a plausible model that satisfies 
both the Herschel and mid-IR data for a CrB has a continuous 
optical depth profile that peaks around 50AU. Because the typi- 
cal expectation from Solar System and protoplanetary disk stud- 
ies is for the surface density to decrease with distanc e from the 
star (e.g. lWeidenschillin3ll977l ; lAndrews et alj|2009h . this pro- 
file can be interpreted as a depletion of material inside 50AU. 
Such a depletion is expected in standard models of collisional 
evolution, where the disk decay rate is faster at smaller semi - 
major axes (e.g. Dominik & Decin 2003; Wyatt et al. 2007a). 
This 50AU turn-over distance is similar to r a ii gn shown in Fig- 
ure^ so it is therefore plausible that the depletion inside here is 
due to increased velocities imposed by binary perturbations. 

Considering this binary stirring picture in more detail, ver- 
tical (inclination) stirring is probably more important than radial 
(eccentricity) stirring in circumbinary disks, particularly when 
the disk lies at large distances relative to the binary separation. 
As the binary mass ratio decreases (i.e. the stars' masses be- 
come more similar), the eccentricities imposed on an exterior 
particle decrease. For example , using the expression derived by 
iMoriwaki & Nakagawal J2004I) . the "forced" eccentricity of par- 
ticles at 50AU around a CrB is about 0.001, at least an order of 
magnitude lower than the eccentricities thought to exist in ob- 

6 Once catastrophic collisions are occurring, the level of stirring is also 
important. An increase in stirring does not have a significant effect on 
the collision rate between particles of similar sizes because higher ve- 
locities are accompanied by a lower volume density of particles (e.g. 
the scale height and/or radial disk extent increases with the velocities). 
However, particles are typically destroyed in a collision with a much 
smaller particle, whose size is set by the collision velocity. Thus, higher 
levels of stirring mean that smaller particles are capa ble of destroying 
objects. In standard collisional size distributions (e.g. Dohnanyj l 19691) 
the number of particles, and therefore the number of potential destruc- 
tive impactors, increases strongly with decreasing size, so the disk de- 
pletion rate due to collisions increases with increased stirring. 
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serve d debris disks (e.g. IKenvon & Bromi ev 2004; Krivov et alj 

2008) . In contrast, for inclination stirring to be effective the ini- 
tial disk plane would only need to be a few degrees different to 
the binary plane (though some non-zero eccentricity is needed 
to ensure crossing orbits). For example, with an initial misalign- 
ment of 1°, collision velocities are stirred to about 0.02 times 
the Keplerian velocity (150m s _1 at 50AU), so objects with dis- 
persal thresholds less than roughly 10 4 J kg -1 (<10km) will be 
disrupted and dispersed in collisions with similar sized objects. 
Therefore, the disk decays from the inside-out as secular pertur- 
bations increase the relative velocities in particle collisions by 
aligning the disk, and do so on a shorter timescale closer to the 
central star(s). Assuming that the initial misalignment was large 
enough, collisions will be destructive inside some radius that 
will b e close to r a ii gn (but not exactly at r a ii gn , iMustill & W vatt 

2009) . Outside this radius, collisions are unaffected by the pres- 
ence of the binary and therefore not destructive, providing a 
possible explanation for the drop in optical depth in our model 
beyond 50AU. 

While the a CrB disk is consistent with being stirred by 
the binary, we also consider two alternative stirring mecha- 
nisms. The first is "pre- stirring" in which objects are assumed 
to have been stirred at some very early time by an unspeci- 
fied mechanism that does not necessarily still operate (e.g. a 
stellar flyby or the result of gas disk dispersal). A potential is- 
sue with such a scenario is that collisional damping may re- 
duce the velocities sufficiently that collisions are not catas- 
trophic for a significant fraction of the stellar main-sequenc e 
lifetime fe.g. IKenvon & Bromlevlbool iGoldreich et al.ll2004h . 
Whether collisional damping is important depends on the rela- 
tive sizes of the objects and their destructive impactors, which in 
turn depends on the relative velocities. If pre- stirred objects can 
reach sufficiently high eccentricities (>0.1, though the value de- 
pends on object strength and stellocentric dis tance), then the 
disk c an remain stirred for the stellar lifetime JShannon & Wul 
1201 ll) . Such a disk is depleted radially from the inside-out, 
again because the depletion rate is a strong function of ra- 
dius. The type of structure that is expected is therefore a ra- 
dially increasing optical depth profile, which turns over and 
decreases w here collisions have yet to deplete the disk signif- 
icantly (e.g . Kennedv & Wvatt 2010). Using equation (6) from 
iKennedv & Wvattl (12010) we find for the best fit planetesimal 
properties from ( Wvatt et al. 2007b) that the a CrB disk would 
be depleted within 50AU if the disk was stirred from an arbitrar- 
ily early time (e.g. pre-stirred) and was about 5-15 times less 
massive than the s olid component of the minimum-mass So- 
lar nebula (MMSN, Weidenschilring l 1977L when scaled linearly 
with the binary mass)Q 

The second mechanism is "self-stirring", where random 
velocities are initially slow enough that collisions result in ac- 
cretion and growth. Once the largest objects reach roughly 
Pluto- size, they increase the velocities of smaller planetesimals 
and initiate a collisi onal cascade, and the p roduction of visi- 
ble levels of dust (iKenvon & Bromlevl 120041) . Again, this pro- 



7 We are not necessarily suggesting that the solid mass in the pri- 
mordial protoplanetary disk was similarly depleted relative to the 
MMSN, as the mass could for example have gone into building plan- 
ets. However, given the large observed dispersion in protoplanetary disk 
masses in star forming regi ons, such a depletion is easily possible (e.g. 
lAndrews & Williams 2 0051) 



cess works in a radially inside-out fashion, so a self-stirred disk 
looks similar to a pre-stirred one, with one key difference. Be- 
cause planetesimals have not been stirred outside where Pluto- 
sized objects have formed, collisions do not result in a colli- 
sional cascade and th e disk should show a drop in optical depth 
beyond this distance (IKennedv & W vatt 2010). Using equation 
(9) from IKennedv & Wvattl ( l2010l) . for Pluto-sized objects to 
stir the disk to only 50AU by 350Myr, the disk would have 
to be 2000 times less ma ssive than a scaled MMSN (see also 
IKenvon &"B romlev 2008, from which the Pluto-formation and 
stirring times were derived). A disk with s uch low mass would 
not be visible in a self-stirring scenario (e.g. lKenvon & Bromlevl 
120081) . which appears to disfavour this scenario. However, this 
calculation assumed that the planetesimals started out with lm- 
lkm sizes. If the planetesimals were initially much larger, the 
time to form Pluto- sized objects would also be longer (e.g. 
IKenvon & Bromlevl [20 10). This longer formation time would 
mean that a more massive disk, which would be correspond- 
ingly brighter and therefore more consistent with the observa- 
tions, could form Pluto-sized objects and stir the disk only as 
far as 50AU by 350Myr. 

Therefore all stirring models appear consistent with the ob- 
served peak in surface density at 50AU. However, the model de- 
rived in ^4.21 required a drop in surface density beyond 50AU, 
a feature expected in self-stirred and binary- stirred disks, which 
would appear to disfavour a pre-stirred interpretation, but not 
discern between self and binary- stirring. A caveat on this con- 
clusion is that the disk structure is poorly constrained by the 
low resolution of the observations, and the dust observed with 
Herschel may not trace the parent body locations, particularly 
in the outer regions. For example, the peak at 50AU may sim- 
ply represent the outer edge of a parent body disk that is pre- 
stirred, and the (decreased) emission beyond 50AU could arise 
from small grains originating at 50AU f orced onto larger ec- 
centric orbits by r adiation pressure (e.g. iThebault et al]|2003l ; 
iKrivov et al1l2QQ6h . 

Regardless of the stirring mechanism, we can compare the 
disk structures with those expected if they decayed from some 
arbitrarily large level. In this picture the face-on geometrical op- 
tical depth (with the same assumptions used above) is 



t — 8.9 x lQ-\ 7/3 M; z/6 L-° h- 1 



(1) 



where M* is th e stellar mass in Solar u nits and t is the sys- 
tem age in Myr IKennedv & Wvattll201QL equation 8). The r 7/3 
dependence gives the expected radial profile of a disk whose 
planetesimal properties are the same everywhere. With these as- 
sumptions, for a CrB the expected optical depth at 350Myr is 
1.15 x 10 -8 at 1AU, which is somewhat smaller than the model 
value of6.4xl0 -8 . Given that there is considerable uncertainty 
in the planetesimal properties and the true dust distribution, we 
do not consider this difference a cause for concern. The model 
does not increase as strongly with radius as Equation (Q]), so is 
below the expected level outside a few AU anyway. That the 
disk model has a r 1 7 dependence rather than r 7//3 could indi- 
cate that the planetesimal properties have a radial dependence, 
for example that they become smaller or weaker at larger dis- 
tances (the model parameters are very uncertain however, so 
such a dependence is not required). 

The f3 Tri data are consistent with a continuous optical 
depth profile that peaks around 100AU. Given that the secular 
precession time depends strongly on semi-major axis, and that 
the stellar age is uncertain, this distance is not significantly in- 
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side r a ii gn and the disk is therefore plausibly binary stirred. The 
expected optical depth at 730Myr is 3.8 x 1CT 9 at 1 AU, which is 
smaller than the model value of 2. Ox 10 -8 . Again, the discrep- 
ancy is not particularly large given the model assumptions and 
uncertainty. Using the same equations as above, the P Tri disk 
would be depleted out to 100AU by 730Myr for a disk 1-5 times 
less massive than a scaled MMSN if it were stirred from the out- 
set. Stirring by Pluto-formation out to this distance only requires 
a disk 400 times less massive than a scaled MMSN (but again 
could be more massive if planetesimals are larger). Thus, the 
disk could be depleted out to 100AU by collisional evolution, 
but the stirring mechanism is unclear. Unlike a CrB, which has 
a mid-IR detection, there is no evidence for warm emission in 
the P Tri system. Such a detection could break the degeneracy 
in our models, which cannot tell whether the regions inside 50- 
100AU are devoid of debris (e.g. due to dynamical clearing by 
planets), or simply depleted by collisional evolution. 

This "standard" picture of debris disk stirring and evolu- 
tion is not the only possibility. For example, catastrophic colli- 
sions may be caused due to crossing orbits in a disk where self- 
gravity is important, with the additi onal possibility that suc h 
disks may appear non-axisymmetric Jjalali & Tremainell2012h . 
It is also possible that some observed debris disks are not stirred 
to cat astrophic collision velocities at all. iHeng & Tremaind 
12013) show that long-lived "warm" planetesimal disks could 
exist, in which collisions are not typically disruptive. They sug- 
gest that a test for such a scenario is that the disk spectrum 
should look like a blackbody, because the disk particles required 
for such warm disks to survive are large enough that they act 
like blackbodies (i.e. absorb and emit light efficiently). For the 
two systems considered here the disk spectra appear to rule out 
such a scenario because they lie significantly below pure black- 
bodies beyond wavelengths of about 100/im (Fig.O, suggesting 
that the particles are emitting inefficiently at long wavelengths 
and are predominantly smaller than ~lmm. However, we do 
not exclude the possibility that the disks are "warm", because 
the small grains that are observed could have been created in 
erosive and bouncing collisions. Modelling the size distribution 
would make more quantitative predictions to test this possibil- 
ity. 

Though both disks have extended or continuous dust dis- 
tributions as observed by Herschel, the parent bodies may (or 
may not) be more localised. All the models we considered find 
that the dust optical depth decreases with distance in the outer 
regions. Qualitatively, this structure is expected when the parent 
bodies occupy a narrow "birth" ring a nd small grains are placed 
on eccentric and hyperbolic orbits fe.g. lStrubbe & Chiand2006l : 
iMulleret all bold) . However, it is also possible that the ob- 
served extent reflects the underlying parent body distribution, as 
might be argued for a CrB, which has dust detected both near 
and far from the star. Alternatively, the drop in optical depth be- 
yond 50AU required by the continuous a CrB model may be a 
sign that the parent planetesimals lie relatively close to the star 
and that the more distant dust comprises small grains on eccen- 
tric and hyperbolic orbits. 

Whether debris disks are typica lly rings or more extended 
structures is an open question (e.g. iKalas et al. 2006), m part 
because obtaining sufficiently high resolution sub-mm observa- 
tions, those most sensitive to larger grains, is challenging. This 
ambiguity has only been overcome in a few nearby systems, 
where a parent body ring has been resolved at sub-mm wave- 
lengths and is seen to be narrower than the radial extent of small 



grains (IKalas et al.ll2004l : IWilner et al]l2012l : iBolev et allEnl 
lAckeetal.ll2012h . With the development of facilities such as 
the Atacama Large Millimeter Array (ALMA) and the North- 
ern Extended Millimeter Array (NOEMA), the detection and 
resolution of larger parent body populations will become pos- 
sible and will provide insight into the processes that set where 
planetesimals form and reside. 

Looking at the issue of alignment from a wider perspec- 
tive, coplanarity is the expected outcome for debris disks and 
planetary systems emerging from the protopla netary disk phas e 
for small to medium binary separations fe.g. lBate et al.| [2000). 
However, only a few examples where the outcome can be tested 
actually exist. Three systems with transiting circumbinary plan- 
ets, in which the stars are also eclipsing bin aries, show that 
well aligned syste ms are a possible outcome ( Dovle et aLll201ll ; 
IWelsh et ahEoIl) . Further. I Welsh et all <2012l) find that the oc- 
currence rate of aligned circumbinary planets is probably con- 
sistent with the rate for circumstellar planets with similar prop- 
erties, suggesting that alignment is the typical outcome. How- 
ever, they acknowledge significant biases exist, and that further 
work is needed to understand the implications of these discover- 
ies for circumbinary planetary system alignment and frequency. 

In the case of circumbinary debris disks, only four sys- 
tems where the disk and b inary alignment can be tested exist; a 
CrB, Tri, 99 Herculis ([Kennedy et al.ll2012h . and HD 988OC0 
(iBoden et al.ll2005l : lAndrews et alJl2Q10h . Of these, a CrB and 
P Tri are close binaries with periods of several weeks and HD 
98800 has a period of 314 days, while 99 Her has a semi-major 
axis of 16.5AU and a 56 year period. With only these systems 
we cannot yet be sure of what trends will emerge. As hinted by 
the disk-binary alignment of a CrB, P Tri, and HD 98800, and 
the misalignment for 99 Her, it may be that more widely sepa- 
rated systems are more likely to be misaligned. It will also be 
interesting to test whether stirring by secular perturbations from 
binaries is important for disk evolution. This hypothesis could 
be tested by comparing disk sizes (and structure where possible) 
with the radii at which secular perturbations can have reached 
within the system lifetime for a larger sample. 



7 SUMMARY 

We have presented resolved images of debris disks around the 
nearby close binary systems a CrB and P Tri. These sys- 
tems are relatively unusual among binaries because their or- 
bital configurations are relatively well known, allowing a test 
for (mis)alignment between the disk and binary orbital planes. 
In both cases we find that the disk and binary are most likely 
aligned. Though secular perturbations can align systems over 
time, the bulk of the resolved disk structure in these systems is 
too distant to be affected. Therefore, the alignment is most likely 
primordial, suggesting that the binary + protoplanetary disk sys- 
tem from which the debris disk emerged was also aligned. These 
initial conditions are consistent with expectations of alignment 
in protoplanetary disk+binary sys tems where the b inary has a 
separation less than about 100AU teate et alfcoOOl) . 

Secular perturbations from the binary could provide the 



8 A weak accr etion signature h as recently been detected for ~10Myr 
old HD 98800 lYang et al.ll2012l) . suggesting that it lies somewhere be- 
tween the protoplanetary and debris disk phases. 
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stirring mechanism in circumbinary disks, and both a CrB and 
f3 Tri are consistent with such a picture. However, they are also 
consistent with other stirring mechanisms. While binary stirring 
may happen, it cannot be the only mechanism because debris 
disks are observed wi th a similar frequen cy in both single and 
multiple star systems (Trill ing et al.ll2007h . 

These two systems bring the number in which debris disk- 
binary alignment can be tested to four. Three of these (a CrB, 
/3 Tri, and HD 98800) have orbital periods of several weeks to 
a year and appear to be aligned, while 99 Her has a period of 
50 years and is strongly misaligned. It is too early to draw con- 
clusions about typical outcomes in such systems, but the results 
so far suggest that misalignment cannot be extremely rare, and 
may be preferred in systems with wider binary separations. 
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